Compressed sensing for multidimensional electronic spectroscopy experiments 
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Compressed sensing is a processing method that significantly reduces the number of measure- 
ments needed to accurately resolve signals in many fields of science and engineering. We develop 
a two-dimensional (2D) variant of compressed sensing for multidimensional electronic spectroscopy 
and apply it to experimental data. For the model system of atomic rubidium vapor, we find that 
compressed sensing provides significantly better resolution of 2D spectra than a conventional discrete 
Fourier transform from the same experimental data. We believe that by combining powerful resolu- 
tion with ease of use, compressed sensing can be a powerful tool for the analysis and interpretation 
of ultrafast spectroscopy data. 



Multidimensional spectroscopy [ll-Ql is an important 
tool for studying ultrafast dynamical processes in com- 
plex molecular systems. For instance, it can be used to 
analyze vibrational energy transfer at liquid/air inter- 
faces on picosecond timescales Q or exciton dynamics in 
natural light harvesting systems at hundreds of femtosec- 
onds [6H8|- It can also be applied for efHcient detection 
and identification of molecules, which is one of the crucial 
challenges in chemistry, biology, and medicine with im- 
portant applications, for example, in molecular sensing, 
chemical separation, and DNA analysis. Frequently, in 
these nonlinear optical techniques, the data collected in 
the time domain is Fourier transformed to the frequency 
domain. The crucial issue, then, for obtaining high fre- 
quency resolution in measured spectra is the long sam- 
pling time required. Here, we demonstrate that a method 
known as compressed sensing (CS) can be applied as a 
very efficient alternative to the Fourier transform to ob- 
tain high-resolution multidimensional spectra. 

Compressed sensing is a state-of-the-art signal process- 
ing method which has recently become popular through- 
out the physical and biological sciences. The method is 
founded on the concept of sparsity. When a signal is 
known to be sparse in a certain basis (i.e. most of the 
coefficients are negligibly small), this additional knowl- 
edge can be used to dramatically reduce the number of 
measurements required to reconstruct the signal [l^l . 
This method has been applied to many areas of research, 
ranging from magnetic resonance ima ging ll| to super- 



resolved imaging of sing le molecules [12| and quantum 
process tomography [13[. Earlier, some of us showed that 
CS can also be used to significantly reduce the compu- 
tational cost of atomistic simulations [l^ . In that work, 
the application of CS was particularized for molecular dy- 
namics and real-time time-dependent density functional 
theory simulations for obtaining linear spectra (vibra- 
tional, optical absorption, and circular dichroism) [Tsj . 
However, CS can be applied to other types of simu- 
lations and experimental techniques. Multidimensional 
non-linear spectroscopy is then an interesting and rel- 



evant candidate to explore the possibilities of CS. Re- 
cently, CS has been pursued to reconstruct one of the 
dimensions in ultrafast 2D NMR data llGll and to simu- 
late 2D spectra using random sampling [17| . 

In this letter, we present the first application of CS to 
experimental ultrafast 2D optical spectroscopy. We ap- 
ply the method to an atomic system - Rubidium vapor, 
which is frequently used as a test model for multidimen- 
sional spectroscopy techniques [l8l - [20j . We find that CS 
presents several additional advantages beyond a simple 
speed-up. Many multidimensional experiments are in- 
herently limited in the amount of time-domain data that 
may be collected, either due to measurement constraints 
or to more fundamental limitations such as the timescale 
of the dynamics one wishes to explore (which may be 
very short due to decoherence and other processes). We 
believe CS can extend the range of spectroscopically- 
observable dynamics by providing higher frequency reso- 
lution even given limited time domain data. 

Another benefit of CS is that it is quite easy for ex- 
perimentalists to integrate. As discussed below, our ap- 
proach simply replaces the 2D discrete Fourier transform 
(FT) by a new 2D compressed sensing scheme. Since we 
recast 2D CS as a series of one-dimensional problems, 
we are able to harness the power of parallel computing 
so that the entire signal processing does not take signifi- 
cantly longer than 2D FT. It is our hope that CS's easy 
portability and short processing time will help it become 
a method of choice among experimental spectroscopists. 

The rest of this letter is structured as follows. We 
begin by presenting the compressed sensing method and 
outlining its application to the resolution of 2D ultra- 
fast optical spectra. Wc next apply the method to a 
model experimental system, namely gas-phase Rubidium 
atoms, and show how CS may be used to better resolve 
the spectral features. We then discuss the experimental 
methods and numerical implementation. Finally, we offer 
conclusions and a future outlook. We begin by describ- 
ing how the CS method can be a ppli ed to 2D ultrafast 
optical spectroscopy; see refs. [21-23| for more detailed 
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FIG. 1. Schematic illustration of the sequence of four pulses 
used in 2D optical spectroscopy. 



information about CS. In particular, we will focus on 
four-wave mixing experiments, similar to those used to 
study coherent energy transfer in light-harvesting com- 
plexes, quantum dots, and other systems of physical and 
biological interest. Typically, these experiments involve 
irradiating a sample with four optical pulses and varying 
the time gaps between the pulses (Fig. [Ij . Defining the 
time gaps as r (the coherence time between pulses 1 and 
2), T (the population time between pulses 2 and 3), and 
t (the waiting time between pulses 3 and 4), the signal 
measured in the time domain is a function S{t, T, t). The 
standard approach in four-wave mixing is to perform a 
2D discrete Fourier transform in r and t to obtain the 
spectrum S{ojT,T,ujt). This spectrum is typically plot- 
ted in ojr — oJt space for different values of T, and the 
dynamics of the peaks then give information about the 
dynamics of the underlying system. 

The CS method can be used to replace a ID discrete 
FT in converting time-resolved data into the frequency 
domain [^, , obtaining a significantly better frequency 
resolution than the conventional discrete FT. Our basic 
approach here, then, is to reduce the 2D CS computation 
to a series of ID CS calculations. A major advantage of 
this approach (as compared to an inherently 2D method) 
is that the ID computations can be performed in parallel, 
and each is individually inexpensive. To obtain the 2D 
spectrum for a particular population time T, our proce- 
dure involves two steps: 

Step 1: For every value of r in the data set, perform a 
ID CS in t to convert the signal to the oJt domain. By col- 
lecting all the results for different t together, we obtain 
the half-transformed signal S{t, T,LUt) which is sparse in 
LOf . This step requires a total of N-r ID CS calculations 
(where Nr is the number of r data points) . Each calcula- 
tion takes a short time and all of them can be performed 
in parallel. 

Step 2: Now for every value of oJt in the half- 
transformed signal, perform a ID CS in r to convert the 
signal to the uir domain. This step requires N^^^ ID CS 
calculations which also take a short time and may be 
performed in parallel. Collecting all the results together 
yields the final sparse spectrum S{uJr,T,ujt)- 

It should be noted that the roles of r and t may be 
interchanged to yield similar final results. In either case, 
if the spectrum is sparse in the LJr — plane, then this 
2D compressed sensing procedure yields better-resolved 
peaks in the frequency domains with less time-domain 
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FIG. 2. Energy level diagram for atomic *^Rb vapor. 



data than the conventional 2D discrete FT, as we will 
illustrate. 

For completeness, however, we first provide a very brief 
review of the ID CS calculations, as these underlie the 
2D method discussed above (see ref. [ij] for a more de- 
tailed derivation). In each ID CS calculation, we want 
to obtain a vector g of values {gi , g2, ■ ■ ■ ■, gN^ } ^'t iV^^ 
equidistant frequencies Uj = Aujj, from the known vec- 
tor h set of time-resolved values {hi,h2, ■ ■ ■ ,hNt} given 
at Nt equidistant times tj — Atj. (The scheme can be 
generalized for non-uniform sampling.) Our objective is 
to obtain sensible results with Nt as small as possible, so 
we are interested in the case N^j > Nt- In principle, the 
gk set can be directly obtained using the discrete FT, 
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However, if we expect that many of the Fourier coeffi- 
cients are negligible, we can use CS to extract more pre- 
cise results from the same signal. This is done by the 
recasting the Fourier coefficient calculation as a linear 
equation. Since this is an underdetermined problem (as 
> Nt), the sparsity condition means that we should 
select the solution that has the larger number of zero co- 
efficients. In practice this solution can be obtained b 
solving the basis-pursuit de-noising (BPDN) problem 
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subject to \Fg — hl^ < rj 



(2) 



where F is the iV^^ x Nt inverse Fourier matrix with entries 
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and rj represents a level of noise that we assume is present 
in the signal. In all calculations which follow, 77 < 10"'*. 

To illustrate the utility of CS in two-dimensional op- 
tical spectroscopy, we consider phase-modulation 2D flu- 
orescence spectroscopy (PM-2DFS) data collected from 
atomic * ""Rb vapor [20| ■ The ^ ''Rb system may be consid- 
ered as a quantum three-level system with ground state 
5 ^S'1/2, first excited state 5 and second excited 
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state 5 ^^3/2, as illustrated in Fig. [2l The four light 
pulses are produced by a titanium saphire laser [with 
full width at half maximum (FWHM) - 42 fs] which 
is on resonance with the two electronic D line transi- 
tions of ^''Rb: 5 ^<S'i/2 5 (with transition fre- 
quency 2.370 rad/fs) and 5 ^51/2 5 ^^3/2 (with tran- 
sition frequency 2.414 rad/fs). We considered ®^Rb va- 
por to be an ideal candidate for our 2D CS method as 
its 2D electronic spectrum is expected to be sparse in 
frequency space, with diagonal and cross peaks corre- 
sponding to the transitions just mentioned. The natu- 
ral lifetime of electronic excitations in *^Rb is ^ 25 ns, 
which should correspond to very narrow lines (D/ 36.1 
and D// 38.1 rad/ms). Full and extensive details of 
the PM-2DFS experimental method used to collect the 
data are given in refs. 
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The experiments yield 
time-resolved fluorescence-detected "sum" and "differ- 
ence" signals [SsumiTjT^t) and S'diff (t, T, t)] which are 
analogous, respectively, to the nonrephasing and rephas- 
ing third-order polarizations collected in more traditional 
four- wave mixing experiments, minus the effects of non- 
resonant interactions. 

The main results of this letter are given in Figs. [3] and 
131 which compare the performance of the conventional 
discrete Fourier transform (top) and our 2D compressed 
sensing method (bottom) in resolving the peaks corre- 
sponding to transitions in the atomic ^''Rb vapor. Fig. [3] 
shows the "sum" (nonrephasing) signal Ssum{'-^T,T,ujt) 
while Fig. |4] shows the "difference" (rephasing) signal 
'S'diff(wr, i^t) for the population time T ~ 140 fs. Sim- 
ilar results were obtained for other population times. 
Note that the exact same set of time-resolved data was 
used for the discrete FT spectra as for the CS spectra. 

As can be seen from the figures, CS produces peaks 
that are far better resolved in frequency space than those 
obtained by the discrete FT. In fact, the peaks obtained 
by CS are at least 5 to 10 times less wide in each di- 
mension than those obtained via the discrete FT, con- 
sistent with the results we previously obtained in one- 
dimensional spectra [14| |. In other words, to achieve a 



comparable resolution as CS with the discrete FT, one 
would need to collect far more time-resolved data by con- 
tinuing to measure the signal S{t, T, t) for longer time 
delays r and t. Given that CS obtains narrow peaks 
with only limited time-resolved data, we expect this sig- 
nal processing method to be particularly useful in resolv- 
ing closely-spaced peaks in their 2D spectra. 

One immediate question arising from the comparison 
of the FT and CS spectra in Figs. [3] and |4] is that of 
peak shape. As for long sampling times the CS result 
matches the FT spectra, the question is how capable is 
CS of resolving those peak features for short times. The 
issue of peak shape certainly deserves a more extensive 
investigation in the future, most likely by applying our 
2D CS method to more complex experimental systems 
with intricate features that are broadened by internal 



structure or an environment. 

In conclusion, we have demonstrated the first applica- 
tion of compressed sensing to two-dimensional electronic 
spectroscopy experiments. Focusing on electronic tran- 
sitions in an atomic ®^Rb vapor model system, we have 
shown that the 2D CS method that we have developed 
provides much finer resolution of the peaks in a 2D spec- 
trum as compared to the standard discrete FT. As a re- 
sult, we expect that 2D CS will substantially reduce the 
experimental effort needed to obtain well-resolved spec- 
tra by decreasing the amount of time- resolved data which 
must be collected. Furthermore, we hope that 2D CS will 
increase the range of spectral features resolvable in ultra- 
fast experiments, particularly closely-spaced peaks. 

An important question is the degree to which CS ac- 
curately reproduces peak shapes in more complex sys- 
tems with internal structure and an environment; this 
question definitely merits a further joint theoretical and 
experimental effort. We hope and expect that CS will be- 
come more widely investigated and employed in the 2D 
experimental ultrafast community once its easy portabil- 
ity, competitive speed, and strong resolving power be- 
come widely known. 



METHODS 

The experimental methods used to collect the time- 
resolved data for the atomic ^^Rb vapor are de- 
tailed in ref. Time-resolved "sum" and "differ- 
ence" signals were collected at all points on a two- 
dimensional grid consisting of 51 equally-spaced coher- 
ence times T and 50 equally-spaced waiting times t (with 
At ^ At = 26.687 fs) for a series of 5 population times 
r= 140, 175, 210, 245, 280 fs. 

The full 2D CS spectrum was obtained by performing 
a series of ID CS calculations as described in the text. 
For each ID CS calculation, the optimization problem in 
eq. © was solved using the SPGLl algorithm developed 
by van den Berg and Friedlander [25| . To avoid numerical 
stability issues we work with a normalized BPDN prob- 
lem, where the prefactor 2 Aw /tt of the F matrix, eq. ([3]) , 
is left out and h is normalized. The missing factors are 
included in g after the solution is found. 

For the CS calculations, we use a frequency grid con- 
sisting of 1000 evenly-spaced points between — tt/At and 
tt/At for ojr and 1000 evenly-spaced points between 
—ir/At and n/At for wt (the same grid spacings are also 
used for the discrete FT). 
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FIG. 3. Comparison between discrete Fourier transform (top) and compressed sensing (bottom) for the "sum" (nonrephasing) 
2D optical spectra of *^Rb vapor for population time T — 140 fs. CS yields narrower, better-resolved peaks than the discrete 
FT. 
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FIG. 4. Comparison between discrete Fourier transform (top) and compressed sensing (bottom) for the "difference" (rephasing) 
2D optical spectra of *^Rb vapor for population time T = 140 fs. CS yields narrower, better-resolved peaks than the discrete 
FT. 
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